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Abstract. We use the Villain approximation to show that the Gross- Witten model, in the 
weak- and strong-coupling limits, is related to the unitary matrix model that describes U{N) 
Chern-Simons theory on S^. The weak-coupling limit corresponds to the q —>■ 1 limit of the 
Chern-Simons theory while the strong-coupling regime is related to the q ^ limit. In the 
latter case, there is a logarithmic relationship between the respective coupling constants. We 
also show how the Chern-Simons matrix model arises by considering two-dimensional Yang- 
04 ■ Mills theory with the Villain action. This leads to a U{1)^ theory which is the Abelianization 

^ — ^ I of 2d Yang-Mills theory with the heat-kernel lattice action. In addition, we show that the 

' character expansion of the Villain lattice action gives the q deformation of the heat kernel as 

, it appears in g-deformed 2d Yang-Mills theory. We also study the relationship between the 

^j' unitary and Hermitian Chern-Simons matrix models and the rotation of the integration contour 

, in the corresponding integrals. 

m 

1. Introduction 

In the 1970s, Wilson introduced and studied lattice versions of Yang-Mills (YM) theory, 
■ in which the dynamical variables are elements of the gauge group, defined on the links that 

connect the adjacent lattice points [1]. The choice of the lattice action in Yang-Mills theory 
, was an important aspect of the problems considered in the early development of lattice gauge 

' theory |l]-[9]. In particular, in the late 1970s and early 1980s, a number of papers studied this 

CN ' possibility, leading to the consideration of several alternatives to the Wilson action [1], like the 



(N 



heat-kernel action [21 [3lll] or the Manton action [5j . One of the motivations was the need for a 
proper understanding of the transition from strong to weak coupling in lattice gauge theories. 
I A well-known result in the study of non-Abelian two-dimensional Yang-Mills theory with the 

Wilson action is the third-order phase transition, found by Gross and Witten [lOj and Wadia 
|11| . in a one-plaquette model, described by a unitary matrix model. This result has turned out 
to be of relevance in many current problems in theoretical physics, like the study of Hagedorn 
^ . and deconfining transitions in weakly coupled Yang-Mills theory [12j. The Gross- Witten model 

I has been also recently discussed in the study of type OB and OA fermionic string theories [1^ 

and in relation with other solvable models, like the Kontsevich model |14j . 

We will first show that the unitary matrix model of U{N) Chern-Simons theory on [151 
[T6| [T7] is intimately related to the Gross- Witten model when one considers it together with 
the Villain approximation of the XY model [181 119j. We introduce the Villain approximation, 
together with the two relevant matrix models, in the next section. We will show that both the 
weak-coupling and strong-coupling regimes of the Gross- Witten model can be described, using 
the Villain approximation, by analytically continued U {N) Chern-Simons on S^. Recall that the 
Chern-Simons action is given by 



(1.1) Scs{A) = ^ [ TiiAAdA + lAAAAA), 

4vr Jm -3 

where A is the connection, a 1-form valued on the corresponding Lie algebra, and G Z is the 
level. The (^-parameter is defined in terms of the level k hy q = exp (27ri/(k + N)). 
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We will see that the weak-coupling limit corresponds to the g — )• 1 limit of the Chern-Simons 
theory, whereas the strong-coupling limit of the Gross- Witten model corresponds to the opposite 
limit, q ^ 0. In the matrix model formulation, the q parameter is treated as real and written in 
terms of a coupling constant gs as q = e"^" [21] ■ The above characterization of Chern-Simons 
theory as being analytically continued precisely refers to this treatment of g as a real parameter. 

Indeed, actual computations with the matrix model are carried out with q real, using for 
example the associated g-orthogonal polynomials [21], and the identification gs = 2iTi/{k + N) 
at the end, allows us to make contact with the well-known expressions for the Chern-Simons 
observables |20j . See for example |21J, where the simple case of the U{N) Chern-Simons partition 
function on is computed with the Stieltjes-Wigert polynomials with a q parameter, q = e~^*. 

The correspondence between the Gross- Witten model and the Chern-Simons matrix model is 
of a rather different nature in the two opposite limits q — )• and g — )• 1. In the weak-coupling limit 
g\y[ — ?• 0, as shown already in [22], the coupling constants are related by ^ym ~ ^^s, whereas 
in the strong-coupling limit we will have that g'g = 21n ((7ym) '^^^ equivalently, q = I/q'ym' 
shall see in Sec. 2 in detail. 

We will end Sec. 2 by exploring some consequences of these relationships between the models. 
In particular, in Sec. 2.1 we show that the Gross- Witten model at weak coupling is a Gaussian 
matrix model whose free energy has an expansion that can be interpreted in terms of closed 
strings. 

The relationship between the Gross- Witten model and the Chern-Simons matrix model, based 
on the application of the Villain approximation to the Wilson action, indicates that the direct 
consideration of the Abelian Villain action in lattice two-dimensional Yang- Mills theory, should 
describe Chern-Simons theory. Two-dimensional Yang-Mills theory was also studied with the 
Manton action [231 124j and the heat-kernel action |4]. We shall see that, indeed, the straightfor- 
ward generalization of the Abelian U{1) lattice action, which is just a theta function [181 El E] 



to the non- Abelian case, in the setting of U{N) two-dimensional Yang-Mills theory, directly 
gives U{N) Chern-Simons theory on S^. This straightforward extension of the Villain lattice 
action to the non- Abelian case was explored by Onofri, shortly after the study of the heat-kernel 
case |3], in a less well-known work |9]. The description of pure Chern-Simons theory by such a 
model has not hitherto been realized. 

It is well known that the Villain model arises in the Kogut-Susskind Hamiltonian lattice gauge 
theory [25] in the Abelian case, which leads to a direct correspondence with the planar Heisenberg 
(or XY) model ^61 12]. The non- Abelian case leads to the heat kernel [21 H], and we shall see 
that the Chern-Simons matrix model follows from Abelianization of the heat-kernel propagator 
in the context of two-dimensional Yang-Mills theory. This is the content of Sec. 3 and, in 
particular, we show in Sec. 3.1 that this Abelian projection is equivalent to a g deformation of 
2d Yang-Mills, in consistency with the known relationship between Chern-Simons theory and 
a q deformation of 2d Yang-Mills theory [2^. Recall that Chern-Simons theory is known to 
be explained in terms of an Abelian two-dimensional Yang-Mills theory, as was shown at the 
level of the path integral, first in the case of manifolds of the type x T,h, where E^, denotes 
a Riemann surface of genus h |27j and, more recently, for Seifert fibrations over S/j |28] . which 
contains the case, the one studied in this paper at the level of the matrix model. 



(1.2) 




l=—oo 
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To conclude, we study in the Appendix the precise relationship between the unitary and the 
Hermitian versions of the Chern-Simons matrix model focusing also in the rotation of the con- 
tours of integration. 



2. Gross- WiTTEN model and the Villain approximation 

The approximation devised by Villain in 1975 in the study of the two-dimensional XY model 
[18j is based on the simple observation that the term exp (/3 cos 9) that appears in the 2d XY 
model can be well approximated for large /3 by a periodic Gaussian with minima in the same 
locations and with the same curvature. That is 

oo 

(2.1) exp (/3cos0) ~ e'^ ^ ^-\P(e-2.nf p ^ oo. 

n=—OD 

But the l.h.s. term is of course also the weight function of the matrix model description of the 
one-plaquette model of Yang- Mills theory based on the Wilson action (namely, the Gross- Witten 
model [IQj). The r.h.s is a theta function and then the Villain approximation applied to the 
Gross- Witten model leads to the relationship with a unitary matrix model with a theta function 
as weight function. 

Precisely, the unitary matrix model that describes U{N) Chern-Simons theory on 5*^ [IS [Ml 
[T7] is given by 

(2.2) 4f^(5^) = / n?0(«''^i«) nh''^--'''i'' 

•^0 j=l k<l 

where the weight function of the matrix model is a Jacobi third theta function 

oo 

(2.3) co{e) = Q{e''^^q)= q^'/V"''. 

n=— oo 

We show now how this model follows from the Gross- Witten model |10^ \TT\ , by using the Villain 
approximation [19]. Aspects of the relationship between the two models, especially in the weak- 
coupling regime, have already been studied in [22], using orthogonal polynomials. The results 
of the seminal works \18\ [T9] also allow us to extend the relationship between the Gross- Witten 
and the Chern-Simons models to the strong-coupling regime. 

Recall that the Gross- Witten model is a unitary one-matrix model which arises as the one- 
plaquette reduction of the combinatorial quantization of Yang-Mills theory. In two dimensions 
the reduction is exact and described by the partition function [TU] 



(2.4) ^^^(/3) := / dU exp Tv{U + U^)) 

JU(N) / 

r27T N 



,2 

sir 

i<j 

where dU denotes the bi-invariant Haar measure for integration over the unitary group U{N). 
The P parameter is usually written in terms of the gauge coupling constant as /? = [lO] 
and hence, the strong coupling limit is given by /3 — )• and the weak coupling limit by /3 — )• cxd. 



-'^The unitary matrix model H2.2|) also describes Chern-Simons theory if the weight function is 0^^( — e |g) 
[16| . This possibility has also been noticed in [T7]. The nonuniqueness description of the Chern-Simons matrix 
models is described in |5T]. See the Appendix for its precise relationship with the Hermitian matrix model 
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The planar or XY model is characterized by a coupling between nearest-neighbor spins which 
has the same analytical form as the potential in the Gross- Witten model [18\ [T9] 

(2.5) V,3 = -/3 [1 - cos {9 - 9')] . 

The coefficients of the Fourier expansion of this potential are given by 

(2.6) e^(^) = /, (/3) , 

where Iniz) is the modified Bessel function of order n. The partition function of the Gross- 
Witten model can also be written as the determinant of a Toeplitz matrix with ()2.6p as entries 
of the matrix |29j 

(2.7) ZnW)= det [/.-,(2/3)] . 

On the other hand, for U{N) Chern-Simons theory on the corresponding Toeplitz determinant 

is [lell 



(2.8) zg^)(53)= det k_j(g)] , with aj=g^'/2. 

l<i,j<A'^ 

In the weak-coupling limit /3 — )• oo, the Fourier coefficient is 

(2.9) hm e^(^) = e'^'/^/? ^ ^-s^g^/^ 

In this limit, the Toeplitz determinant (j2.7p has the Gaussian coefficients (j2.9p as entries, and 
hence it coincides with (|2.8p . Of course, (^ym ~^ ^ ™ (|2.9p and this implies — )• on the 
Chern-Simons theory side as well. 

The prescription in |18j . namely ()2.ip . is valid for both the opposite /? — t- and /3 — )■ oo limits, 
using always periodic Gaussians, but including a renormalization scale Ry (/3) and a rescaled 
inverse temperature fiy = f {f3) . This implies that a correspondence between the Gross- Witten 
model and the Chern-Simons matrix model holds for both the weak-coupling and the strong- 
coupling regimes of the model. As seen above with the Toeplitz determinant representation of 
the matrix model, and also in [22] using orthogonal polynomials, the weak-coupling limit follows 
in a straightforward way. This result, as shown in |19j . also follows by considering decimation 
[30j . It is explained in |19] that, after a few iterations of the Kadanoff-Migdal decimation 
procedure, any interaction function at reasonably low temperatures generates a new interaction 
of the Villain type 

oo 
m=—OD 

In addition, the critical properties of the two models could be identical by conveniently choosing 
/3v = / (/3) a function of f3 [19j. In particular, a comparison of the interaction form for weak 
and strong coupling showed equivalence if 

/(/3) = /3 for /3^oo 

(2.10) /(/3) = [21n(2//3)]-i for (3^0. 

Notice that to compare with the theta function as it appears in the matrix model (|2.2p . we have 
to take into account that 

oo ^ oo 

(2.11) g-/3(0+27rn)2 ^ ^ ^ Q-nVm^ine^ 



^Such a determinant was already considered in [Oj. See Sec. 3. 
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The r.h.s. of (j2.11|) is the series expansion of the theta function in (|2.2|) . Thus, the decimation 
of the weight function of the Gross- Witten model coincides with the ViUain approximation and 
it leads to the weight function of the unitary Chern-Simons matrix model (j2.2p . In the weak- 
coupling limit /3 — 7- oo, we obtain that gs = 1//3 and hence that gs = g\y^l1. This also follows 
from (E^D. 



The strong-coupling regime /3 — )• is specially interesting. Notice that taking /3 = {g\}^ 
oo) directly leads to the circular ensemble |31j 

277 //). a \ N 



= 0) = ^ ^ n ( V") n • 



In [22], the relationship between this model and Chern-Simons theory was studied. The result 
(j2.10p leads to a refined understanding of this limit /3 — 0. In particular, the relationship 
between the coupling constants is now logarithmic gs = 21n(2//3) = 21n (sym) hence, the 
coupling constant of the Gross- Witten model is related to the q parameter of Chern-Simons 
theory q = l/gy^. To summarize, including the prefactors 



(2.12) e^™^^ « _^e(ei^.|e-V/3), with5s=9W2forgYM^0, 

cose ^ e(e '^^1 e-i/'^v), with g^ = 2log{g^y,) for ^ym ^ oo. 

It is also worth mentioning that, already in the original paper on the XY model [18], an approx- 
imation valid for both limits was also given 



(2.13) e^"""'^ ^ RviP) J2 ^' 



■l3vie-2Tvrnf/2 



with 



(2.14) Rvi^) = /o(/3)V2^ 

6-/^^2 = /i(/3)//o(/3), 

where Ii {(3) and /o(/3) are, as in (|2.6|) . Bessel functions. These are the first two Fourier coef- 
ficients of the expansion of the l.h.s. of ()2.13p . The values for the prefactor Ry {(3) and the 
coupling constant fiy are then found by imposing the first two Fourier coefficients of the two 
expressions in (|2.13p to coincide. The limits /3 — )■ and /? — >• oo of (|2.14p coincide with the 
previous results. It was shown in [18] that this approximation is rather good, even for values 
close to the critical temperature. This approximation has been studied in further detail in [32], 
where it was found that the convergence is better for the strongly coupled regime /3 — )• 0. In 
addition, the result ()2.13p can be extended, with good convergence in both limits, to the case 
where the original function is u) (9) = /3cos6 + 7 cos (26) [32j. This suggests that not only the 
Gross- Witten model, but more complex multicritical unitary matrix models [33], can also be 
expressed in terms of the unitary Chern-Simons matrix model, only with a more sophisticated 
expression for the coupling constant and the prefactor. 



2.1. Gaussian behavior and closed string interpretation at weak coupling. The re- 
lationship between the Gross- Witten model in the weak-coupling limit and the semiclassical 
limit of the Chern-Simons matrix model, indicates that the Gross- Witten model in this regime 
should be related to closed topological strings |34j . The reason is that the semiclassical limit 
{k — )• 00) of the U{N) Chern-Simons on free energy also coincides with the nonperturbative 
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part of the total free energy. The Chern-Simons free energy can be suitably expressed in terms 
of nonperturbative and perturbative contributions 

(2.15) Fes = log Zcs = Fnp + i^p . 

The nonperturbative contribution F^p is the logarithm of the measure factor in the path integral, 
which is not captured by Feynman diagrams, and it gives the exact Chern-Simons partition 
function in the semiclassical limit A: — )• oo |34l eq. (2.8)]. It has the explicit expression 

This free energy ()2.16p , given by a Hermitian Gaussian matrix model, has an expansion which can 
be interpreted in terms of closed topological string theory on the resolved conifold geometry [M] . 
See pi] for equivalent string theory and gauge theory interpretations. 

Due to the correspondence between the Chern-Simons and the Gross- Witten matrix models, 
the latter should have (j2.16|) as free energy in the weak-coupling limit /3 — t- cxd. Consider the 
expression for the free energy of the Gross- Witten model for finite and ^ym — ^ [35j 

2_/y _/Y _/Y^ 2 N—i 

Fgw^^-^ ln(27r) - — In(^) + In j! + ©(^^-m)- 
9ym ^ ^ Sym j=i 

This coincides with the free energy of a Hermitian Gaussian matrix model |31] 



(2.17) Zg (/3) = e^^ / e- M/^ " 



This relationship between the Gross- Witten and the Gaussian matrix model agrees with |36) . 
Notice that the term e'^^ actually corresponds to the e^ term in the Villain approximation (j2.12p 
and implies that one has to consider the Gross- Witten model with the potential , instead 

of the one in S22 



3. The Villain lattice action and Abelian/^-deformed 2d Yang-Mills theory 

We begin by discussing the heat-kernel action, which was introduced in lattice gauge theory, 
at least in part, as an alternative to the Wilson action that also provided a natural extension 
to the non-Abelian case 121 [3l H] of the Villain approximation of the two-dimensional XY model 
[15] . which had been crucially used in the study of C/(l) lattice gauge theories [SI [7]. 

The setting is quantum Yang-Mills theory with gauge group U{N) on an oriented closed 
Riemann surface S/j of genus h and unit area form dfi |37j . The action is 

(3.1) 5ym = -t^ / d/iTrF^ , 

45YM JEfe 

where Qs plays the role of the coupling constant, F is the field strength of a matrix gauge connec- 
tion, and Tr is the trace in the fundamental representation of U{N). A lattice regularization of 
the gauge theory relies on a triangulation of the two-dimensional manifold S with group matrices 
situated along the edges [38]. The path integral is then approximated by the finite-dimensional 
unitary matrix integral 

(3.2) Zm = J Yl dUe J] Zp[Up] , 

edges I plaqucttcs P 

where dUi denotes Haar measure on SU{N) and the holonomy Up = II^gp ordered 
product of group matrices along the links of a given plaquette. The local factor Zp [Up] is a 



This form of the potential corresponds exactly to the Wilson lattice action. 
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suitable gauge invariant lattice weight that converges in the continuum limit to the Boltzmann 
weight for the Yang- Mills action ()3.ip . This expression for Zm leads to the matrix models 
presented in the previous section after the suitable choice of lattice action. 

Let us discuss now the use of the heat kernel for the lattice weight Zp [Up]. In this case, 
the lattice action has many interesting features [4j and is the usual choice in two-dimensional 
Yang- Mills theory [S^. It leads to the well-known group theory expansion of the partition 
function 

(3.3) = ^ (dim A)^^^'' exp ( - gyu C2W) , 

A 

where the sum runs through all isomorphism classes A of irreducible representations of the 
SU{N) gauge group, dim A is the dimension of the representation A, and C2(A) is the quadratic 
Casimir invariant of A. This expression for the partition function is a particular case of the 
propagator [H [37] 

exp (-5hk m = (I| exp f ^^^-mA^ \U) = KiU, %), 



which can be written as [1] 

(3.4) K{U, ^) = J^dimAxA (U) exp {-g^yM>^)), 

A 

where the sum runs over all irreducible unitary representations of the gauge group, xx (U) is 
the character of such a representation, dim A = xx (I) its dimension and C2(A) the Casimir of 
the representation. Writing (|3.4p in terms of the elements of the Young tableaux that labels the 
representation A leads to a discrete matrix model representation [37] . 

3.0.1. q deformation. In recent years, the relationship between two-dimensional Yang-Mills 
theory and Chern-Simons theory on Seifert manifolds has been understood in further detail 
p6l HOl E] (see also [28l HH US])- Seifert manifolds M{h,p) are nontrivial circle bundles (of 
monopole degree p) over two-dimensional surfaces of genus h. The simplest case, the trivial 
fibration, M(/i, 0) = T,h x was studied in detail in \20\ 127]. 

Chern-Simons theory on Seifert manifolds has been the subject of much interest in the study of 
topological strings |44] and has a direct relationship with a q deformation of the two-dimensional 
Yang-Mills theory discussed above, when the manifold is a sphere, Sq = 5"^. In particular, the 
partition function of g-deformed 2d YM on a closed Riemann surface of genus h is given by [26] 

(3.5) Z^,,(S,) = j;(dim,(A))^-^Ni^^W, 

A 

where dim^ (A) is the q deformation of the dimensions of sin representations, i.e. the quantum 
dimensions dimg(A) [IS], p is a positive integer parameter and, as usual, C2(A) is the Casimir 
of the representation A. This is related to the partition function of Chern-Simons theory on 
a circle fibration (with Chern class p) over S/j, which is a Seifert space. In the case = S'^ 
and p > 1, the Seifert manifold is the lens space S^/Zip. If p = 1, then the connection is with 
Chern-Simons theory on S^, the case studied here. 



3.1. Abelianization and q deformation. The propagator (|3.4p can be alternatively written 
in terms of the elements of U{N), and then one obtains a unitary matrix model expression. 
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When written in terms of the invariant angles of the gauge group, it is given by |1] 
(3.6) 



exp(-5HK(^i,...^7v)) =AA n 



{0= 



-oo i<j 



gj + 27r jlj-lj) 
2 sin [9i - 9j + I-k {U - Ij] 



■ exp 



N 



9ym 



where N stands for some normahzation. This is a rather complex model, and therefore the heat- 
kernel case, in contrast to the two cases studied in the previous section, is not studied with a 
unitary matrix model but rather with a discrete matrix model that follows from (|3.3p . We show 
now that a simplification of ()3.6p leads to the Chern-Simons matrix model. The exponential 
part in the r.h.s. is the Abelian U{1) Villain action, which is just a theta function (|2.ip 



(3.7) 



exp(-Sv(0))= « 



/= — OO 



and the straightforward generalization to U (N) gives a propagator 



(3.8) 



1 ^ 

— Y{ej + 27Ti,] 



9ym 



{/}=_oo 
N 

JJexp(-S'v(6'i)) = exp(-S'viiiain(6'i,---6'Af)) 



1=1 



that defines an action which is just the direct product of the Abelian Villain action ()3.7p . The 
corresponding one-plaquette model is given by a U{N) matrix integral, which is the unitary 
integration of the propagator. The partition function is then given by the matrix model 



(3.9) 



2n N 



da 



nil E 



-(6»i+27m)^ 



n 

k<l 



which, recalling the identity (|2.1ip . is the unitary matrix model description of U{N) Chern- 
Simons theory on ()2.2p discussed so far and previously studied in [151 US [IZl- Hence, 
Zj\[ = Z^^^^ {^^)^ by identifying gs = Qym/'^- Thus, the Abelianization of the heat kernel 
()3.6p . given by the straightforward extension of the U{\) Abelian Villain action to the U{1)^ 
case p.Sp . leads to U{N) Chern-Simons theory on S^. 

At the level of the matrix model, this straightforward generalization of the Villain action 
was already considered in an interesting paper by Onofri [9], that analyzed this problem and 
presented a detailed study of the matrix model (j3.9p . It is remarkable that the results in that 
paper describe pure Chern-Simons theory on S^, with U{N) and SU{N), a fact that has not 
been pointed out so far. Some expressions in [9| are manifestly Chern-Simons observables, like 
the partition function of U{N) Chern-Simons theory on S^. 

We have seen the Abelianization of the heat kernel, using its expression in terms of the 
invariant angles of the gauge group ()3.6p . Let us now look at the relationship between the choice 
of the Villain lattice action and the q deformation of the heat-kernel 2d Yang-Mills theory, using 
instead the character expansion ()3.4p . 

Notice that the g-deformed 2d Yang-Mills (j3.5p is just given by the usual expression for 
two-dimensional Yang-Mills theory based on the heat-kernel action and defined on a compact 
manifold of genus h, but with q dimensions instead of ordinary dimensions, dim(A) — )• dimq(A). 
We will now show that choosing (j3.8p with the Villain action (jl.2p as a lattice action is actually 
equivalent to g-deformed 2d Yang-Mills. Since choosing the Villain action is tantamount to an 
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Abelian projection of the two-dimensional Yang-Mills theory with the heat-kernel lattice action, 
this implies that the Abelianization and the q deformation of the latter lead to the same theory. 

An expression of the type (j3.5p is valid for 2d Yang-Mills on a compact manifold of genus 
h and it can be constructed from the central heat kernel, which is also the propagator on a 
cylinder [37] 

(3.10) K{^,U,U') = J2xx iU)xx {U') exp (-^^^.CslA) ), 

A 

with holonomies U and U' on the two disks of the cylinder and where xx (U) is the character 
of the irreducible representation A, as in (|3.4p above. Recall that xx {U = 1) = dim(A). Hence, 
ii U =[/' = !, ()3.10p is the expression for the partition function of two-dimensional Yang-Mills 
theory on S"^. If only U' = I, then it describes two-dimensional Yang-Mills on a disk 

(3.11) K{^, U,I) = Y^ XX (U) dim A exp {-g^y^,C2{X)). 

X 

This is the heat kernel and it is also the amplitude of a plaquette that leads, by gluing, to 
()3.6p [37 1 . Since (|3.1ip is also the fundamental solution of the diffusion equation, then, as is 
well-known, two-dimensional Yang-Mills theory based on the heat kernel can be understood as 
a diffusion process on the gauge group manifold. 

The relationship between the Chern-Simons model constructed with the Villain action and 
the q deformation of two-dimensional Yang-Mills theory can be readily seen by considering the 
character expansion of the U{N) Villain action (|3.8p given in [9] 

e-'S'viiiain(c/) 1 jv 2 T-r /a'^'""'^ 1\ 

(3-12) "-^ = E X(n.....,n..)(^)g^^-"'n( 

mi>m2>...>mjv=0 j>i 

where the sum is over the integers {rrii} with i = 1, and Z is the partition function of the 
unitary matrix model, which has the explicit form 

2 \ f I. 2 
9ym\ TT ^^ ^N~kY 



^ k=i 



2 



This character expansion of the Villain action suggests that such an action leads to the usual 
propagator of two-dimensional Yang-Mills theory based on the heat kernel ()3.1ip but with q 
dimensions instead of dimensions, since the last term in (j3.12p gives an explicit expression for 
quantum dimensions. An elementary manipulation of the product shows this explicitly 

j>i ^ ^ ^ j>i Iq 

where A denotes the unitary irreducible representation of the gauge group, characterized by a 
partition whose Young tableaux has columns with {mj} boxes. The explicit expression for the 
quantum dimensions is given, as in [26], by 



TT L^-i ~ "^7 + J ~ 'Ja 
dimg(Aj =11 j— — Tj with the (7- number pjq = 

j>i ~ — 



1 1 



Putting all together, we see that the Casimir term is manifest and appears exactly as in [26] 

N N N 

C2(A) = E + (^ - 2« + l)mi = Yl "'^("'^ -2i + l) + N] 



=1 
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and hence we obtain 



e 



■^VillainCC^) 



(3.13) 



Z 



X 



Thus, the r.h.s of (j3.12p is the disk amphtude for g-deformed two-dimensional Yang-Mihs theory 
with p = 1, which, by the same procedure discussed above, and explained in detail in [37J, leads 
to the expression for the partition function ()3.5p . This identity is actually a particular case of 
the Kostant identity, which gives a character expansion of theta functions of a lattice [36] . This 
specific form of the identity was then rediscovered later on, in [9], in the context of the Villain 
lattice action, and it also appears much later in [17], where, working from the r.h.s. of ()3.13|) a 
theta function expression was found, which is, as we have seen here, the Villain lattice action. 
In this paper, we have also shown, by using the equivalent formulation of the lattice action in 
terms of the invariant angles ()3.6p . that it follows from taking only the Abelian part of the heat 
kernel. 

The fact that the Abelianization and the q deformation described here are equivalent is 
qualitatively consistent with the fact, explained for example in [IS], that a q deformation of a 
Lie group G is not a group and lacks its symmetry, whereas the maximal torus T of G remains 
an ordinary symmetry group after the symmetry breaking inherent in the transformation of a 
Lie group into a quantum group. 

To conclude, let us mention that inspection of the analogous expression for SU (N) in [9j 
shows that dim^ A can also be written as the character x\ (Tq), where Tq £ SL (N, C) is given by 
diag q^~^, ...,9^"^, • Hence, the propagator is now of the type K{^, U, Tq) instead 

of (|3.1ip . which suggests that diffusion does not take place in the whole gauge group. Indeed, 
the explicit form of the matrix model indicates diffusion on the maximal torus U{1) x ... x U{1) 
of the gauge group U (N). This is in agreement with a previous result that related Chern-Simons 
theory on 5'^ with Brownian motion on the Weyl chamber of the gauge group |49j . 



We have seen how the unitary matrix model that describes U (N) Chern-Simons theory on 
arises from studying two-dimensional Yang-Mills theory with the Villain lattice action and 
we have compared it both with the Wilson and the heat-kernel lattice action cases. 

Regarding the former, we have seen that the Gross- Witten model is related to the Chern- 
Simons matrix model both in the weak-coupling and the strong-coupling regimes. As we have 
seen in Sec. 2.1., one of the implications is that the Gross- Witten model, which describes 
two-dimensional Yang-Mills theory on M^, coincides in the weak-coupling limit with the non- 
perturbative part of Chern-Simons theory on S^, and consequently has the same string theory 
interpretation [M]. In both cases, the free energy is given by a Hermitian Gaussian matrix 
model. In spite of the apparent simplicity of such a matrix model, it is actually relevant in 
the study of subsectors of AA=4 super symmetric gauge theory and their relationship with two- 
dimensional Yang-Mills theory (see [50], for example). It is possible that taking into account 
Wilson loops in our discussion would lead to a relationship with that line of research. 

Notice also that the general approximation of the Gross- Witten weight (|2.13p implies that a 
small modification of the Gross- Witten model potential leads to the Chern-Simons matrix model. 
The interest of this result lies in a possible connection between the Chern-Simons matrix model 
and the unitary matrix models that appear in the study of phase transitions of weakly coupled 
gauge theories [T2] . 

The heat-kernel has a character expansion which is the basis of the study of 2d Yang-Mills 
theory. However, it also can be expressed in the invariant angles of the gauge group (the unitary 
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group in our case), as pointed out in [4j. The unitary matrix model that fohows from this 
representation is not used in the heat-kernel case, due to its complexity. However, we also 
have seen that an Abelian projection of the heat-kernel lattice action leads to a U{1)^ lattice 
action, which is the Villain lattice action. After unitary integration of the resulting propagator, 
the corresponding matrix model is now the U{N) Chern-Simons matrix model for S^. On the 
other hand, since the character expansion of the Villain lattice action gives the g-deformed 
2d Yang-Mills propagator (Kostant-Onofri identity), we see that the Abelianization of (|c{.6|) 
coincides with the q deformation of ()3.1ip . given by ()3.13p . In addition, it shows how the q 
propagator directly leads to the unitary Chern-Simons matrix model, instead of the (equivalent) 
Chern-Simons Hermitian matrix model. 

Precisely, and to conclude, in the Appendix, the relationship between the unitary and the 
Hermitian versions of the Chern-Simons matrix model is studied in detail, focusing also in the 
rotation of the contours of integration. 

Acknowledgments. Thanks to Sergio Iguri and Fokko van de Bult for comments and corre- 
spondence. The work of MT has been supported by the project "Probabilistic approach to finite 
and infinite dimensional dynamical systems" (PTDC/MAT/104173/2008) at the Universidade 
de Lisboa. 
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Appendix A. Unitary vs Hermitian and integration contours in matrix models 

In this paper we have focused on the unitary matrix model that describes Chern-Simons 
theory on . It was first considered in [15] although we have seen that it was already studied in 
detail in ^ . In contrast to the Hermitian matrix model, no additional work has been done using 
the unitary model, with the exception of its recent appearance in the study of matrix models in 
Donaldson-Thomas theory [IJl fT6| . The Hermitian matrix model, in the case reads [51] 

(A.l) Zg^) {S') = r e- Yl ( 2 sinh " V ;Q ^ . 

Both models, the unitary (j2.2p and the Hermitian (jA.ip . are solved with orthogonal polynomials 
that have the same orthogonality properties, and hence the corresponding observables of the 
matrix model, like the partition function for example, coincide. These polynomials are the 
Stieltjes-Wigert polynomials in the Hermitian case [21] and the Rogers-Szego polynomials in 
the case of the unit circle [52^ I22j . The connection between these two systems of orthogonal 
polynomials allows to explain the relationship between both matrix models, as was shown in 
|52j . However, it is interesting to have a more immediate relationship between the models. Let 
us write the model (|A.ip in its trigonometric version: 

2 N 

(A,2) is') = r e-i^." n (2-^) n 

j<k ^ ^ j=i 

We will now show its relationship with (jA.ip . The first step is to relate (|A.2p with the unitary 
model (j2.2p . by transforming the Vandermonde determinant 



(A.3) ]J Asin^^^^y = []|e*«^-e^'^^P, 

j<k ^ ^ j<k 

and the weight function, making the range of integration compact and using also the identity 
(I2TT]1 



j<k j=i (271") 2 Jo 27r 

r27T N 



JV 

2 



e 



iMfe 1 2 



Hence, we see that the trigonometric matrix model ()A.2p is equivalent to the unitary Chern- 
Simons matrix model ()2.2p . It is also related to (lA.ip by the change of variables x = iu 

(A.5) e- ^ n f 2 sinh f ^i^I^) ) ' dx, 

However, notice that the integral in the trigonometric model (jA.Sp is actually over the imaginary 
line of the complex plane. Then, to show the equivalence of (jA.ip and the unitary model (j2.2p . 
the r.h.s of (|A.5p has to be equal to (|A.2p . What we get now is (jA.2p but with the opposite 
sign in the weight and a different integration contour. If we were able to rotate the contour to 
the real axis, then the sign in the exponential of the weight can be corrected just by complex 
conjugation of (|A.5p . recalling as a purely imaginary quantity. Hence, one needs to rotate 
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the contour to the real line and due to the results above mentioned, we know that this has to 
be the case. 

This rotation of the contour is also directly related to the fact that the actual derivation of the 
matrix model for Seifert manifolds leads to integral expressions with contours on the complex 
plane [53l [51], and it is just assumed that they can be rotated into the real axis |5ij, leading 
then to expressions such as (jA.ip . 

It is then worthwhile to examine the rotation of contours more carefully. For this, we use 
certain particular cases of multivariate hyperbolic hyper geometric integrals, studied in [54j, that 
we can easily identify with the Chern-Simons matrix integrals but with complex integration 
contours. The integrals in [53] are written in terms of the hyperbolic Gamma function. However, 
due to the following property of the hyperbolic Gamma function [5^ 

(A.6) —— = -4sin( — ) sin( — ) , 

it is then immediate to identify such integrals in \5^ with those given by the matrix model 
description of Chern-Simons theory on S^. To see that this is the case, let us first give some 
definitions in [54j and then we quote Proposition 5.3.19 in [53]. 



Definition (hook): A hook W^-^^^p^ is a contour parametrized by W^-^^^^i^) — se^'^^ for s E 
(—00,0] and W^i,<^2(*) = se*'^^ for s € [0,oo). The following quantities are defined as well 
04- = max(arg(a;i), arg(a;2)) and 0_ = min(arg(a;i), arg(a;2)) and, in terms of 4>j^ and (p-, the 
domain A-\- in the complex plane 



-4+ = {4>+-Tr 
After using ()A.6p . Proposition 5.3.19 reads 
Proposition 1: For t G Z>o we have 



+ 0_ - vr. 



(A.7) = ^ e 6*-i-2 ll(^2sin(^^j 

where the contour of integration C is a hook W^^^^j with 0i, (/>2 G A+. 

This result can be extended to t G Z<o by using a property of invariance under complex 
conjugation 



(A.8) JN,t{0Ol,UJ2) = JAf-t(-t^l, -W2). 

To compare with the expression given by the Chern-Simons matrix model, we take oji = L02 = i', 
therefore = (— fjO). We also perform the change of variables 27rxj = Ui and identify 
t = 2{k + N). Then the integral can be written as 
(A.9) 

J^llismhH^^^)e^^^-^ll^ = iV!e-^(fe+Ar)-f e-^^^ {[ (2 sin (^))''^' 

j<k j i — 1 
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or, equivalently [using (jA.Sp ] 
(A.IO) 

/ n 4 sinh2(^^^)e-^ u^Yl^ = Nle-^ (fe+Ar)-f e^^^ {[ (2 sin (-^) ) 

If we move the contour of integration of (jA.lOp to the real axis, we get the expression for the 
Chern-Simons partition function computed in |21j using the Hermitian matrix model. Notice 
that in [53], when restricted to the case S^, the partition function is then given by a single 
complex integral over the line M x e~^, which is exactly (|A.10p for that choice of contour. 

We can use ()A.9p and (jA.lOp to establish the relation between the unitary and Hermitian 
matrix models. If we rotate the contour of ()A.9p to the imaginary axis [that means, taking 
the hook W^^-^{s)], we get the unitary matrix model with q = e~^^, by using the identities 
previously derived 



9 ^ / / • \ \ N-j 

inN^ iV _ 9siV''(iV-l) -p-|- I . I TTJ ^ ^ 

k + N 



(A.ll) = N\e — —{k + N)'-e 12 _Q I 2 

By noting that (jA.lOp is precisely ()A.ip . then 

N 

Jo \ 2^ 7'A \9s J 



j<k 

N-l 

^-i^NiN+Dj^^Ylil-q 



i\N- 

i=l 



This establishes the precise relation between the unitary and Hermitian Chern-Simons matrix 
models, derived directly from contour rotation and is our final result. It is worth mentioning 
that other integrals considered in [54J correspond to the Chern-Simons matrix model on 5^ but 
for other gauge groups . For example. Proposition 5.3.18 in [5lj reads 

Proposition 2: For t £ Z>o we have 

JN,ti^i,^2) ^ -^=-j^ lle-i-2^^dx,llism{^ ^)sin(^^ ^ 



UJl UJ2 

N 



l<r j=l 

J|4sin(-^±- ^Jsinf-^:^^ ^jj_[2sin(^ 



iv -e 

t2 j<k ~ ' ~ ' j=l 

where the contour of integration C is a hook with <j)i,(j)2 G ^+ . This gives an evaluation 

for the partition function of the Chern-Simons matrix model for and the symplectic gauge 
group Sp{2N). 
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